Tunable refraction in a two dimensional quantum metamaterial 
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In this paper we consider a two-dimensional metamaterial comprising an array of qubits (two level quantum 
objects). Here we show that a two-dimensional quantum metamaterial may be controlled, e.g. via the application 
of a magnetic flux, so as to provide controllable refraction of an input signal. Our results are consistent with a 
material that could be quantum birefringent (beam splitter) or not dependent on the application of this control 
parameter. We note that quantum metamaterials as proposed here may be fabricated from a variety of current 
candidate technologies from superconducting qubits to quantum dots. Thus the ideas proposed in this work 
would be readily testable in existing state of the art laboratories. 
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Quantum metamaterials, i.e., artificial optical media, which 
maintain quantum coherence over the signal traversal time, 
hold promise of becoming a testing ground for the inves- 
tigation of the quantum-classical transition, interesting new 
phenomena in wave propagation, and unusual technological 
applications |fl]|2]. With strong analogies existing between 
atomic physics, quantum optics and superconducting systems 
it is natural to seek technologies that span these fields J3). 
Solid state quantum metamaterials is one such class of de- 
vices where such parallels can be leveraged to great utility. 
Indeed, emphasising such a synergy, the implementation of a 
quantum metamaterial in the optical range is feasible HUD. 
In our view, the experimental realisation of the concept is 
likely to be achieved first in the microwave range, as was 
the case with conventional metamaterials [5]. We believe that 
the best candidate system system would comprise supercon- 
ducting qubits J6] |7) playing the role of controllable artificial 
atoms. This view is supported by the ability of superconduct- 
ing flux qubits to play the role of quantum scatterers - a phe- 
nomena that has been both theoretically modelled and experi- 
mentally observed |[8ilT0t 

One dimensional quantum metamaterials can be realised 
readily enough. An example of this is a chain of qubits placed 
in a transmission line 0^8). Unfortunately, ID-devices do 
not allow us to realise more interesting and useful effects such 
as "quantum birefringence" and other phenomena, where the 
change of direction of the signal by an arbitrary angle is im- 
portant (e.g., [4 1). Here one needs to go beyond ID. As with 
any truly quantum circuit, we retain the essential requirement 
that any "proof of principle" realisation of a truly quantum 
metamaterial must maintain global quantum coherence, and it 
is therefore necessary for such a system to contain as few unit 
elements as possible. In this work we are therefore concerned 
with the minimal realisation of a 2D quantum metamaterial. 

The effect of changing signal transmission amplitude 
through a system by tuning the quantum states of the con- 
stituents, which we consider in this paper, is related to two 
classes of phenomena: mesoscopic transport in quantum in- 



terferometers (for a review see, e.g.,Ch. 5 of ifTTI or Ch.3 
of 0) and electromagnetically induced transparency [EIT], 
effect known in optics lfT2Ul4l and recently demonstrated in 
an artificial atom (superconducting qubit) in the microwave 
range [ 10]. The underlying physics of either of these phenom- 
ena is the same: quantum state-dependent interference be- 
tween different quantum trajectories contributing to the prob- 
ability amplitude of the signals transmission through the sys- 
tem. In the case of EIT, in the presence of one ("coupling") 
laser field destructive interference occurs between different 
transitions in an atom under the influence of the probe laser 
field; as a result, the field effectively does not interact with 
the atom. The control is exercised through the coupling laser 
field, which produces the appropriate state of the atom; the 
mechanism works for each atom in the medium indepen- 
dently. 

In case of mesoscopic transport through a quantum interfer- 
ometer, the transmission amplitude is determined by the inter- 
ference of electron wave components propagating through the 
branches of the device. Here the interference pattern occurs in 
real space, and it can be directly affected by the electromag- 
netic field acting on the electrons most spectacularly, through 
Aharonov-Bohm or Aharonov-Casher mechanisms |fT31 . 

The propagation of electromagnetic signals through a pla- 
nar quantum metamaterial, a toy model of which we are con- 
sidering here, is determined by quantum interference in real 
space, like in mesoscopic transport. On the other hand, one 
cannot directly affect the phase of the electromagnetic wave 
at any reasonable field magnitude; instead, like in EIT, we 
control the coupling between, and thus setting the "preferred" 
states of, the artificial atoms, which constitute our medium 
(with the advantage that here in principle quantum state of 
each of them can be controlled individually and at will). 

Systems of coupled spin-^ particles have been widely stud- 
ied in a range of different physical contexts. Recently, there 
has been much interest in the potential use of spin chains or 
networks as buses for quantum state transfer within quantum 
information processing devices lfl6l IT7l . It has been shown 
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FIG. 1: (color online) Set up, each node represents a qubit and each 
edge the qubit qubit coupling coupling. The directed couplings cor- 
respond to Sij ((f)) = ofaj exp[i<jf>] +a~a^ exp[— i<j!>] and the undi- 
rected couplings to Hij (O). 



that the phase shift caused by an applied field (the Aharonov- 
Casher effect lfT31 ) can significantly enhance the maximum 
attainable degree of entanglement in a spin chain, as well 
as improving the transfer of entanglement around a ring of 
spins 1181 . The effects of continuously monitoring the out- 
put nodes of a spin network was investigated in [19] and it 
was demonstrated that high-fidelity quantum state transfer is 
possible. We note the emphasis of [19| differs from ours as 
they were concerned with communicating a particular state 
while we are interested in characterising a particular mate- 
rial. Indeed, we do not ask questions about the actual state of 
the system as a whole which will mostly be some entangled 
state of the qubits; we are simply interested in an output signal 
comprising local measurements on two edges of the material. 

Our simulated quantum metamaterial comprises an array of 
interacting qubits as depicted in Fig.[T] In this work we con- 
sider the reduced model for qubits coupled to electromagnetic 
modes where the field degrees of freedom have been elimi- 
nated producing effective couplings between the qubits simi- 
lar to [20]. The Hamiltonian for this system is then 
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where erf is the Pauli matrix for the z direction and we have 
chosen jUjj = g for the connected qubits as illustrated through 
the use of edges in Fig. [T] and zero otherwise, i.e. A = 

{(1, 2), (2, 5), (5, 4), (4, 1), (2, 3), (5, 6), (5, 8), (4, 7), (7, 8), 
(8, 9), (9, 6), (6, 3)}. The coupling operator Sjj(0jj) takes 
the form: 
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where af — \ (erf ± icrf ). In order to model the effect of 
(for example) a magnetic field on this qubit metamaterial we 
have chosen the control parameter, cf>ij = <fi for the directed 
couplings (1,2), (2,5), (5,4), (4,1) and zero otherwise as 



indicated in Fig. [T] The initial state of the system was cho- 
sen to be | — 1)^ where |— l) i is the eigenstate of erf with 
eigenvalue —1. We note that for our candidate realisation of 
superconducting qubits the phase <fi may be switched at high 
frequencies, e.g. of the order of a few GHz, ET1 . Conse- 
quently, rapid control of a quantum metamaterial such as the 
one proposed here should be possible. As it may be of utility 
when considering classical analogues of this system we also 
note, that it is possible to write the Hamiltonian in the follow- 
ing form: 
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which may be considered as a twisted XY model. 

Our motivation for this particular quantum metamaterial is 
a quantum circuit realisation fabricated using and array of su- 
perconducting qubits. This justifies our chosen Hamiltonian 
as it is of the same form as those currently used to model such 
circuits. It is, however, interesting to note that it would be pos- 
sible to construct an equivalent system of interacting fermions 
(this can be achieved via a Jordan-Wigner or similar trans- 
form 1 22]). The results presented here in terms of the states 
of a qubit can, therefore, also be viewed in terms of the mo- 
tion of fermions. Such a process might, for example, take the 
form of electrons propagating in a lattice of quantum dots in 
the presence of a magnetic flux which is precisely the case of 
mesoscopic transport mentioned previously. 

In order to characterise the properties of our proposed quan- 
tum metamaterial we need to generate some kind of excita- 
tion flow through the material. Hence, we will need an input, 
which we arbitrarily have coupled to qubit 1. Moreover, we 
also need to measure the output of this flow on the outer edges 
of the material, that is, qubits 7, 8, 9, 6 and 3 (reading left to 
right around the opposite edges to qubit 1). We have chosen 
to do this through the introduction of quantum jumps based 
measurement devices (as the outcome of a measurement is 
clear, either a jump is measured or it is not, see for exam- 
ple ||23l ). In effect, we measure in and create excitations on 
node 1 in Fig. [T] and measure out and destroy excitations on 
the opposing edges (nodes 7, 8, 9, 6 and 3). 

The quantum jumps model is an unravelling of the 
master equation corresponding to the irreversible emis- 
sion/absorption of absorbed/emitted excitations over very 
short time scales. The equation for this unravelling takes the 
form of stochastic Ito increment equation for the state vector 
according to 
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where B = {1,3,6,7,8,9} and dNj is a Poissonian noise 
process satisfying dNj dN^ — SjkdNj, dNjdt = and 
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FIG. 2: (color online) We represent along the abscissa each of the 
qubits to which the output measurement devices are coupled, i.e. 
qubits; 7, 8, 9, 6 & 3 (as indicated by the schematic beneath the 
plot). On the ordinate axis we show the phase <j> for the coupling 
term H(</>) of Eq. {2j. In the z direction (indicated through the col- 
ormap) we show the total number of counts in each detector over ten 
different solutions to Eq. {3}. Here we have used the following mate- 
rial and detector parameters; /xy = | (for the connected qubits) and 
Tin = 7out = 1 (we note that results of a similar form are found for 
IMij = -|). 

dN] = (l\L 3 \ dt. That is, jumps occur randomly at a rate 
that is determined by i^L^Lj^. Here the Lindblad associated 

with the input node is L\ = ji n &t and with the output nodes 
is L; = Toutc, 7 for i £ {3, 6, 7, 8, 9} and zero otherwise. 

In order to take account of the statistical nature of unravel- 
ling's of the master equation we have, in the results that fol- 
low, summed the counts measured over ten trajectories. For 



each trajectory we integrated the counts measured over a rea- 
sonable duration of time (0 < t < 1000). These two factors 
together are sufficient to average out most of the statistical 
fluctuations that arise when using quantum trajectories meth- 
ods (we are confident that the slight asymmetries that remain 
in our results are due to the remaining fluctuations). 

On inspection of the results presented in Fig.[2]we see sev- 
eral points of interest. From the coupling term in Eq. |2]) 
we first observe that we would expect our results to be sym- 
metric about <fi — as is indeed the case. If, together with 
this coupling term, we also examine the loop between nodes 
l->2->5->4^1 of Fig. [T] Here the phases are applied to 
impose <fi with a hop along the arrow and —0 in the opposite 
direction. We see that reversing all the arrows is equivalent 
to changing the sign of cf> and equivalent to a reflection of the 
lattice about the axis going through sites 1,5 and 9. Reflecting 
this symmetry we see that the results of Fig. [2] are invariant 
on changing <fi to —<fi and simultaneously swapping labels be- 
tween 3 and 7 and between 6 and 8. 

For certain values of the controlling parameter (phase) 4> we 
see behaviour that is consistent with that of a quantum bire- 
fringent material acting as a beam splitter. We see this from 
4> = ±7r to around ±§7r (where the effect is strongest) where 
we see that the vast majority of counts are measured on qubits 
7 and 3. This pattern is again repeated around <f> = ±j7r. 
Although not conclusive, the possibility that the material is 
acting as a quantum birefringent beam splitter is supported by 
the calculation of the second order correlation coefficient 

(2) = (AN 7 AN 3 ) 
9 (N 7 )(Ns) 

between the 7th and 3rd qubit at = 3ir/4 which we 
have determined to be approximately 0.92 which, by analogy 
with l24l . is indicative of non-classical scattering by the mate- 
rial (here, Ni is the record, as a time series, of photons counted 
at each detector). In contrast, at <j> — we see that counts are 
measured in a bell shape centred around qubit 9. Here it could 
be argued that this behaviour corresponds to propagation of a 
signal as a beam across the diagonal of the quantum metama- 
terial. 

Finally we note that at (f> = ±\k and approximately (f> = 
±§7r the distribution of counts is approximately flat over all 
the output qubits. At these points the material is somewhat 
equally opaque across all the detectors with a quantum ana- 
logue of electromagnetically induced transparency occurring 
for the other possible values of <f>. 

In conclusion, we have presented a model of a two- 
dimensional quantum metamaterial whose behaviour is tune- 
able. For different values of the control phase our results are 
consistent with (i) a quantum birefringent beam splitter (ii) re- 
duced transmission plain wave propagation and (iii) a beam 
directed along the leading diagonal. This system contains 
many degrees of freedom and the analysis we have presented 
here only just begins to scratch the surface of what we be- 
lieve quantum metamaterials will be capable of. Some exam- 
ple applications, that are beyond the scope of this work, might 
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be demonstrating possible violations of Bell's inequalities for 
the material in its "birefringent" state (e.g. at <f> = ±37r/4) 
or further modifying the materials properties by applying the 
control flux in a manner that breaks the symmetry we assume 
in this paper. 
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